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I.  INTRODUCTION 


A simplification  to  the  missile  vulnerability  code  (MV) 1 has  been 
derived  which  takes  into  account  the  linear  nature  of  existing  missile 
trajectory  data.  This  simplification  allows  for  approximately  a fourfold 
decrease  in  the  number  of  parameters  required  to  describe  missile 
vulnerability  data.  The  major  use  of  MV  data  is  in  wargaming  programs 
where  computation  time  and  storage  space  are  generally  quite  limited. 
Thus,  a reduction  in  the  number  of  descriptive  parameters  will  benefit 
wargamers  using  the  data. 


II.  DISCUSSION 


The  MV  code  simplification  is  based  on  two  assumptions.  The  first 
assumption  is  that  the  laser  vulnerability  (LV)  code  generated 
(probability  of  damage)  data  may  be  parameterized  by  a ’’slanted  step" 
function.2  That  is,  the  missile  as  a function  of  laser  dwell  time  t 
may  be  represented  as: 


O 

II 

0 

tth>t 

P = 

P 

/ t-tth  \ 

t , < t < t 
th  max 

*D 

DMAX 

\ tmax~  ^'th  / 

PD  = 

PDMAX 

t < t 

max 

[1] 


where  t , = threshold  dwell  time 
th 

P . . v = maximum  P^  achievable 
DMAX  D 

t = maximum  effective  dwell  time 
max 

The  second  assumption  is  that  the  missile  trajectory  data  may  be 
represented  by  a linear  function  of  miss  distance  versus  initial  damage 
range.  Thus  we  have 

x = mR  + b [2] 

where  x = miss  distance 
R = initial  damage  range, 

M.O.  Sahiaky  "Missile  Vulnerability  (MV)  Code  User’s  Manualy"  Ballistic 
Research  Laboratory  Report  No.  1924.  (AD  #B013530L) 

J.T.  Klopcicy  "A  Proposed  Laser  Vulnerability-Wargaming  Interface" , 
Ballistic  Research  Laboratory  Memorandum  Report  No.  2533  (AD#B0074286) 
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and  m,b  are  constants  determined  by  fits  to  missile  trajectory  data, 
where  damage  has  been  initiated  at  various  points  along  the  trajectory. 
The  above  formula  accurately  describes  missile  trajectory  data  obtained 
to  date. 


The  relevant  output  of  a missile  vulnerability  analysis  is  the 
differential  probability  of  impact  as  a function  of  miss  distance, 
dPj/dx.  For  the  damaged  missile  there  is  a one-to-one  correspondence 

between  the  damage  probability,  P , and  the  impact  probability  Pj. 

That  is,  an  incremental  change  in  the  missile  PQ  implies  an  equal  change 

in  the  impact  probability  P^,  or  dP^/dP^  = 1.  Using  this  simple 

relationship,  the  relevant  quantity  dP^/dx  may  be  derived  from  Equations 

1 and  2.  From  Equation  1 we  have,  for  t ,<  t < t , 

’ th~  max’ 


PDMAX 

t -t  , 
max  th 


[3] 


From  Equation  2 we  have 


dx 

dR 


m 


[4] 


Then,  combining  Equations  3 and  4,  we  get 


dP  = dPl/dt  * J. 

dx  - dx/dR  dR/dt 


DMAX 

Tt  -t  , )m 
v max  thy 


[5] 


where  V = dR/dt  is  the  missile  velocity.  Now  for  most  situations,  the 
missile  velocity  may  be  assumed  constant  over  the  time  interval  from 
tth  to  tmax'  ^us>  we  see>  from  Equation  5,  that  dP^/dx  is  constant 

over  this  time  interval.  A plot  of  dPj/dx  versus  x would  yield  a flat 

distribution  with  limits  defined  by  the  values  of  x corresponding  to 

values  of  R at  the  times  t,  and  t . These  values  of  x are  simply 

th  max  r J 


where 


‘min 

= mR  . + b 

mm 

[6a] 

max 

= mR  + b 

max 

[6b] 

‘min 

= R.  . . - v * t 

initial  max 

[7a] 

‘max 

= R.  . . - v * t, 

initial  th 

[7b] 
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Here  R.  . . ..  is  the  initial  range,  R is  the  range  at  which  the  P„ 
initial  6 * max  6 D 

begins  to  accumulate,  and  R . is  the  range  at  which  P_WAV  has  accumulated. 

It  is  assumed  that  the  dP^/dx  distribution  cuts  off  sharply  at  the 

limits  X . and  X . Realistically,  we  expect  the  distribution  to  fall 
mm  max  ' 

off  more  gradually  at  the  end  points.  The  miss  distance  as  calculated 
from  Equation  2 is  not  completely  deterministic,  and  actual  miss  distances 
would  be  distributed  about  the  calculated  value.  In  reference  1 a 
Gaussian  distribution  about  the  calculated  miss  distance  was  assumed, 
and  the  accumulated  probability  of  kill  was  distributed  over  a range  of 
miss  distances.*  Here  we  ignore  the  tails  of  the  distribution  in  order 
to  simplify  the  results.  The  error  made  in  doing  this  is  negligible, 
as  will  be  shown  in  the  next  section. 


III.  EXAMPLE 

A sample  calculation  is  performed  to  illustrate  the  simplicity  of 
the  derived  methodology.  Let  us  assume  that  the  missile  velocity  is  260 
metres/second  and  the  missile  CEP  (circular  error  probable)  is  7.6  metres. 
Assume  the  miss  distance  equation  is  given  by 

X = . 25R  - 100  [8] 

The  following  quantities  define  the  missile  PQ  (see  Equation  1) : 


t = 2.2  sec 

max 

Pn  = .9 

Dmax 

We  first  calculate  the  miss  distance  distribution  for  a dwell  time  of 
1 second.  From  Equation  5 we  have 


dP  (.9)  1 

dx  (2)  .25  260 


.0069  /metre 


From  Equation  7 we  calculate 

R = 1000-260  * .2 
max  n.0 

= 948.  metres 


[9] 


[10] 


R . 
mm 


1000-260  * 1. 
740  metres 


■k 

See  Appendix  B in  Reference  1 


[11] 
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DWELL  TIME  1 * 


8 12  16  2i 

X , CEP 


Damaged  Missile  Impact  Distribution 


10 


0 4 8 12  16  20  24 

X , CEP 


t igure  2.  Damaged  Missile  Impact  Distribution 
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Thus,  from  Equation  8,  we  have 


max 


= . 25(948. )-100 
= 137.  metres 


X . = .25(740.) -100 

min  = 85.  metres 


The  distribution  will  be  plotted  in  CEP  intervals. 
9,  12  and  13  to  CEP  units,  we  get 


[12] 

[13] 

Converting  Equations 


dP, 


max 


mm 


(.0069) (7.6) 

= .052/CEP 

[14] 

= 137. /7. 6 ~ 

18  CEP 

[IS] 

= 85. /7. 6 ~ 

11  CEP 

[16] 

We  note  that  3 quantities,  dPT/dx,  X , and  X . are  sufficient  to 

describe  the  distribution.  By  contrast,  the  original  formulation  of  the 
problem  required  that  dPj/dx  be  specified  for  each  CEP  interval.  Thus 

to  specify  the  distribution  out  to  10  CEP  from  the  origin  required  10 

numbers,  and  information  about  the  distribution  beyond  10  CEP  was  lost. 

The  distribution  is  plotted  in  Figure  1 as  a solid  line.  Also  shown 

in  Figure  1 is  the  same  distribution  with  the  effects  of  dispersion 

included  (dashed  line).  It  can  be  seen  that  the  tails  of  the  distribution 

are  small  and  can  be  neglected.  The  same  calculation  will  now  be 

performed  for  a 2 second  dwell  time.  We  note  that  the  calculations  of 

dP  /dx  and  X are  the  same  as  before.  To  calculate  X . we  nave 
I max  min 


and  thus 


R . 


min 


1000-260  * 2. 
480.  metres 


X . 


min 


. 25(480. )-100 
20  metres 
20  ~ 3 CEP 

7.6 


[17] 


[18] 


The  distribution  is  plotted  in  Figure  2 (solid  line),  and  again  we  note 
that  the  tails  of  the  dispersed  distribution  (dashed  line)  are  negligible. 
Thus  far  we  have  been  concerned  with  the  impact  distribution  for  a 
damaged  missile.  The  damaged  missile  impact  probability  Pj  is  obtained 

by  integrating  dPj/dx  over  the  miss  distance  X.  In  Figure  2,  for  instance, 

the  integrated  probability  under  the  solid  line  equals  0.78.  There  are, 
however,  two  other  sources  of  probability  for  a given  engagement.  One 
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is  the  probability  of  catastrophic  kill,  P^,  and  the  other  is  the  survival 

probability  Pg.  The  probability  of  catastrophic  kill  is  the  probability 

of  defeating  the  missile  before  ground  impact,  i.e.,  a warhead  kill,  while 
the  survival  probability  is  defined  as  the  probability  of  the  missile 
completing  its  mission  without  being  damaged.  The  sum  of  the  three 
probabilities  must  be  unity  to  conserve  total  probability: 

Pj  . Pc  * Ps  = 1 [19] 

The  survival  probability  Pg  is  distributed  according  to  the  usual 

circular  Gaussian  distribution  (see  Reference  1).  Let  us  assume  for 
simplicity  that  the  probability  of  catastrophic  kill  P„  equals  zero. 

We  then  have 


P 


S 


[20] 


= .22 

The  survival  probability  is  plotted  in  Figure  3 along  with  the  impact 
distribution  from  Figure  2. 


IV.  SUMMARY 

A simplification  to  the  MV  code  has  been  derived  for  situations 
where  the  missile  miss  distance  is  a linear  function  of  damage  range. 
This  simplification  provides  for  a more  concise  representation  of  the 
data  and  is  advantageous  for  wargaming  applications.  A sample  Fortran 
program  which  will  perform  the  calculations  described  in  sections  II 
and  III  is  listed  in  the  Appendix. 
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♦MSTEMP 

.MVFAS 

1 

C 

THIS  ROUTINE  WILL  CALCULATE  THE  IMPACT  DISTRIBUTION  ^OR  A i 

2 

c 

MISSILE.  REQUIRED  INPUTS  ARE  PARAMETRISED  LV  DATA  (SLANTED  i 

3 

c 

STEP  FUNCTION) , THE  MISSILE  CEP  AMO  VELOCITY  , ANO  THE  MISS  j 
DISTANCE  AS  A FUNCTION  OF  RANGE#  THE  IMPACT  DISTRIBUTION 

4 

c 

5 

c 

IS  PRESENTLY  CALCULATED  FOR  RANGES  OF  t '2,3,4,  AND  5 KM, 

6 

c 

AMD  FOR  DWELL  TIME?  OF  1,2*3'  AND  4 SEC#  THE  OUTptIT 

7 

c 

VARIABLES  ARE  THE  NEAPEST  (TO  MISSILE  T*RGrT>  C.rp  INTERVAL 

8 

c 

IMPACTED,  THE  FARTHEST  CEP  INTERVAL  IMPACTED'  AMp  THE  TOTAL 

9 

c 

PROBABILITY  0 e KILL,  WHICH  IS  A FLAT  DISTRIBUTION  BETWEEN 

10 

c 

THESE  TWO  POINTS. 

11 

c 

12 

c 

13 

dimension  XAP (3) 

14 

DATA  CEr>/7.6/VEL0C/260./DMDR/.2S/rON/-lO0./ 

15 

INTEGER  INF/IO/OUTF/II/ 

18 

DATA  RMIN/400./ 

17 

c 

18 

c 

19 

c 

RF AD  IN  LV  DATA  ON  FILE  IN? 

20 

75 

RFAD( INF,100,ERR  = 1000,END=?000)  ID, lNSrO, THETA » PrI , XAP , 

21 

*BRY 

»PKPWRT» SIGMA, Tl, T? »P2»CHIS0 

22 

100 

FORMAT(A2»I4,?F4,0»3F5.1»F4.2»lPER»2»n»F4,l»3X,3FS.2,!DEt‘)o2)‘ 

23 

c 

24 

c 

INITIALIZE 

25 

ICEP1  = 0 

26 

ICEP2  = 0 

27 

PROBK  = 0.0 

28 

IF ( P2 ,LE» 0 , ) GO  TO  25 

29 

c 

30 

c 

LOOP  OVER  INITIAL  RANGE 

31 

DO  1 I = 5, 1,-1 

32 

RANGI  = FLOAT ( I ) *1000 , 

33 

c 

LOOP  OVER  DWELL  TIME  (TD) 

34 

DO  1 J = 1,4 

35 

TP  = PLOAT(J) 

36 

IF(TD.LE.Tl)  GO  TO  25 

37 

c 

COMPUTE  INITIAL  DAMAGE  RANGE 

38 

RT  = RANGI  - VELOC+T1 

39 

c 

IS  T2  GREATER  THAN  DWELL  TIME 

40 

TFIN  s T2 

41 

IF(T2,GT.TD>  TFTN  = Tp 

42 

c 

COMPUTE  FINAL  DAMAGE  RANGr 

43 

RF  = RANGI  - VELOC*TFim 

44 

c 

COMPUTE  ACCUMULATED  PK 

45 

PPOBK  = (TFIN  - Tl)*P?/(T2-Tn 

46 

c 

DETERMINE  cep  interval  limits 

47 

Ip (Rl .LE.RMIN)  GO  TO  ?5 

4R 

XWIS2  = DMDR*RI  + CON  i 

49 

ICEP2  = XMIS2/CEP  + *5 

50 

IF  (RF.GT.RMIN)  GO  TO  10 

51 

c 

SPECIAL  CASE,  DWELL  TIME  UP  TO  RMIn  IS  CALCULATED 

5? 

ICEP1  = 1 

53 

TDD  = (RANGI-RMIN)ZVELOC 

54 

PROBK  = <TD0-T1)*P2/(T?-TD 

55 

GO  TO  2S 

56 

10 

CONTINUF  19 

57 

XMIS2  = DMDR*RF  * CON 

58 

ICEPl  = XMIS2/CEP  ♦ .5 

59 

25 

CONTINUE 

60 

C 

WRITE  DATA  FOR  ONE  RAn®E»  DWELL  TIME 

61 

WRITE (OUTF# 200 #ERR=3000»EnD=4000>  ID* INSEQ# THETA ,PHI» PRY# 

62 

« PKPWRT  * S I SMA  » VELOC » RANOI  » CEP » TO# ICEP1  * IC£P2 » PRO0K 

63 

200 

FORMAT (A2# 14 #2F4.0»P4.1»  lp£8.  ?» 0=>F4. 1 #^4, 0 » F5 . 0 # f4. 1 » F? . n , 

64 

$2 1 3 

»F4.2> 

65 

1 

CONTINUE 

66 

GO  TO  75 

67 

C 

68 

C 

EPROR  EXITS 

69 

1000 

WR ITE ( 6» 1001 ) INF 

70 

1001 

FORMAT ( * ERROR  IN  REApING  FILE* #15) 

71 

STOP 

72 

3000 

WRITE(6»3001)  OUTF 

73 

3001 

FORMAT ( • ERROR  IN  WRITING  FILE*»I5> 

74 

STOP 

75 

4000 

WRITE(6»4001) 

76 

4001 

FORMAT ( » END  OF  FILE  REACHED  ON  WRUe  COMMAND# > 

77 

STOP 

78 

C 

79 

C 

NORMAL  EXIT 

80 

2000 

WRITE (6#  2001 ) 

81 

2001 

FORMAT (*  JOB  COMPLETED*) 

82 

STOP 

83 

END 

RT#S  MSTfWP.WFAS 
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